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The paper offers a novel approach to parameter estimation of a single-diode solar cell/panel equivalent 
circuit, based on analysis of either technical characteristics supplied by the manufacturer or user- 
obtained experimental I-V curve. The derived model allows predicting the solar cell/panel output for 
arbitrary environmental conditions. The method combines a solution of an algebraic equations system 
with an optimization algorithm. The main advantage of the proposed method is that in order to obtain 
the required parameters of the model, minimal set of experimental data is required. In order to validate 
the feasibility of the proposed approach, several solar panels of different types from several manufac- 
turers were analyzed. The results demonstrate 0.1-0.5% estimation precision when the Standard 
Operation Conditions data from the manufacturers’ datasheets are employed. 
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1. Introduction 


Modeling of photovoltaic devices and simulating their behavior 
represent significant portion of today's research in the field of 
solar energy. Modeling using an equivalent circuit (lumped para- 
meters model) is of particular interest, since it allows performing 
joint simulation of a photoelectric device with power electronics 
interfaces [1-5]. Such simulations allow optimizing the design of 
solar arrays and power systems and provide valuable assistance 
during Maximum Power Point Tracking algorithms development 
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without costly experiments but using computer models only [6-8]. 
Hence, proper selection of an equivalent circuit topology along 
with suitable methodology of extracting optimal set of equivalent 
circuit parameters, characterizing a specific device, is a very impor- 
tant part of modeling process. 

Dozens of scientific and engineering publications covering a 
variety of methods for extracting the parameters of equivalent 
circuit models were revealed recently. Most researchers propose a 
topology, including both linear and nonlinear components [9]. The 
linear part consists of resistive elements and irradiation dependent 
current source. Typically, two equivalent resistors are included: 
series and shunt/parallel resistances [10-12]. In some cases, the 
shunt resistance is neglected since it is much larger than the series 
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one [13,14]. One or more ideal diodes typically represent the 
nonlinear part of the equivalent circuit [10,12,15,16]. 

Two fundamental approaches are used for extracting the model 
parameters of a solar cell/panel. The first one involves multiple 
optimization techniques, the most prominent example of which is 
the least mean squares approach. There, the desired parameters 
are obtained by minimizing the difference between measured 
data and the modeled I-V curve. The method requires fluency 
in sophisticated algorithms [17-20,22-24] and returns accurate 
values of the desired parameters in case a sufficiently large 
number of reference points of the measured I-V curve is available 
[21,25-27]. 

Solving systems of algebraic (usually transcendental) equations 
is the second fundamental approach to parameter extraction of a 
solar cell/panel equivalent circuit [28,29]. The amount of experi- 
mental data used in these methods must match the number of 
unknown parameters. The equations are established taking into 
account the constraint that the modeled I-V curve should pass 
through three basic points of the measured I-V curve. Additional 
information such as the derivatives of I-V curve at the basic points 
may be used as well [30-34]. In [13], parameter extraction of a 
single-diode model was performed by solving a system of four 
algebraic equations, three of which represent the I-V curve 
passing through a set of the following basic points: short circuit 
current, open circuit voltage and maximum power. The fourth 
equation used the fact that the power curve derivative is zero at 
the maximum power point. 

The main advantage of analytical methods is that in order to 
obtain the required parameters of the equivalent circuit, minimal 
set of experimental data is required [35]. In many cases, the 
required parameters could be extracted directly from the manu- 
facturer's data. The method in [13] returns good correspondence 
between theoretical and empirical data for single crystalline 
photovoltaic cells and panels for relatively high irradiation levels. 
However, neglecting the shunt resistance causes significant diver- 
gence between the theoretical and real I-V curve for low irradia- 
tion levels. Moreover, this assumption is invalid for modeling 
polycrystalline, amorphous silicon, and organic solar cells under 
any irradiation level. 

In the present study, a novel approach is proposed for extracting 
the single-diode equivalent circuit parameters, combining both 
above mentioned fundamental approaches, i.e. minimization of 
divergences between the modeled I-V curve and experimental data 
by fitting and solution of algebraic equations (partially transcen- 
dental) system. The approach is not new and has been recently 
discussed in [51,52]; nevertheless, unlike in the mentioned refer- 
ences, here the diode ideality is used as the fitting parameter since 
its value resides in a relatively narrow, well-defined region. In 
addition, no simplifying assumptions are made during the analysis. 
The proposed method allows reducing the number of simulta- 
neously solved equations to two. The other two parameters can be 
then obtained by direct calculation, while the fifth parameter is 
derived through optimization/fitting. 

Equivalent circuit of the single diode model adopted in the 
current study is shown in Fig. 1. The model includes a parallel 
combination of a photogenerated current source lpn, a diode D, 
described by the well-known single-exponential Shockley equa- 
tion [36], a shunt resistance R,, and a series resistance R;. An 
arbitrary load is connected to the ports of the model and is 
characterized by the current and voltage I, and V, respectively. 

Note that an analytical solution of system in Fig. 1 does not exist in 
general case. Utilizing the Lambert-W function simplifies the numer- 
ical solution [14,37-39]. This function is the inverse of the function 
f(x) =xe*, ie. W(x) =f~'(x), where W(x) is the Lambert-W func- 
tion. The solution process is consequently simplified since the 
Lambert-W function is a built-in function in most of the existing 


Photovoltaic Cell 


Load 


Fig. 1. The “single-diode” equivalent circuit of the photovoltaic cell with a load. 


engineering software packages, e.g. lambert w in Matlab and Product- 
Log in Mathematica. 

A significant contribution of the proposed method is the 
possibility of obtaining estimated values of band gap energy of 
the cell and the temperature coefficient of the photo-generating 
current. The value of the band gap is usually assumed constant; 
nevertheless, its value varies from one material to another and 
strongly depends on temperature [40]. Moreover, since the single- 
diode model only approximates the underlying physical processes, 
some unmodeled behavior exists in any case and is reflected 
by variations in band gap value. The information regarding the 
temperature coefficient of photo-generating current absent from 
datasheets (the temperature coefficient of the short-circuit 
current is usually given). Most of the researchers utilize the given 
coefficient for estimating the photo-generating current for differ- 
ent ambient conditions by assuming that short-circuit current and 
photo-generating current are equal. Here such assumption is not 
used. Instead, the temperature coefficient of photo-generating 
current is analytically calculated. Equations allowing estimating 
the band gap value and the temperature coefficient of photo- 
generating current require information regarding temperature 
coefficients for both open circuit voltage and short circuit current, 
which are possible to obtain either by test measurements or from 
manufacturer specifications, as mentioned. Thus, only the elec- 
tron's charge and the Boltzmann constant are assumed constant in 
the proposed method. 


2. Materials and methods 


The proposed methodology consists of four main steps. The first 
step involves obtaining the I-V curve expressions based on the 
single-diode equivalent circuit, also referred to as “five parameter 
model”. In the second step, substitution of the standard temperature 
conditions (STC) data into expressions derived during the first step 
creates a set of four algebraic equations, still insufficient for 
determining the five parameters of the chosen equivalent circuit. 
Thus, the fifth parameter is determined by minimization of the 
difference (i.e. by fitting) between the corresponding sets of mod- 
eled and experimental I-V curves in the third step. Here, the diode 
ideality factor is chosen as the fitting parameter [37]. The fourth step 
involves a solution of two additional equations, in which thermal 
coefficients of open circuit voltage and short circuit current are 
utilized to determine the band gap energy and temperature coeffi- 
cient of the photo-generating current. 

It is important to note that all the equations are derived for a 
single while the manufacturer's data is usually given for a panel. 
Since a typical solar panel consists of multiple serially connected 
cells, the panel voltage is scaled down by number of cells to obtain 
single cell data [41]. Nevertheless, in case a solar panel consists of 
parallel cell connection, the panel current should be similarly 
scaled. 
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3. Theory/calculation 
3.1. Deriving the set of equivalent circuit equations 


A set of algebraic equations is derived according to Kirchhoff's 
laws, taking into consideration the characteristic of all circuit 
elements. The diode current may be described by Shockley's 
expression [36] as 


Ip =Ip(e2"" — 1), (1) 


where Vr [V] is thermal voltage, Ip [A] is reverse bias saturation 
current and y is diode ideality factor. The thermal voltage is 
expressed as 


Vr= s Vr (=): (2) 
where kp=1.380 x 107” [J K~*] is Boltzmann constant, q= 1.602 
x 10~'° [C] is elementary charge and T [K] is cell/panel tempera- 
ture. The asterisk (*) here and thereafter denotes a value corre- 
sponding to STC, i.e. irradiance of 1 kW m~?, AM 1.5 and cell/panel 
temperature of 25 °C. 

The reverse bias saturation current depends on temperature as 
[42] 


nk T 3 (Ex/q nVz)(T —T*/T) 
Io =Ig pr) ern ; (3) 


where Eg [eV] is the band gap energy. According to [43-45], and 
[13], the photocurrent is given by 


Toh = Kn pln en(lte(T—T*)) (4) 
where e [K~'] is the temperature coefficient of the photocurrent, 
Y [W m7 ?] is the irradiance and kpn is the spectral factor (kp = 1 
for AM 1.5) [13]. 

Applying Kirchhoffs laws to Fig. 1, the set of equivalent circuit 
equations consists of the three following relations, 


Inn —Ip —Isn —I, = 9, (5) 

Vp—V,—I,R; = 0, (6) 
Vp 

lnr= s 7 

sh Ron (7) 


Note that resistors R; and Rsn as well as the diode ideality factor 
yn are further assumed to be temperature independent [13]. 
Combining (1) and (5)-(7) and applying the Lambert-W function, 
the expression of I-V curve may be presented in the two following 
forms, 


Vi = (lo +lpn) ` Rsn— tn (Rs +Rsn)—1 Ver - wir -Rgp © e0 71+ ph) Ron/01° mi 
MANE 


(8) 
and 


(lo +Ipn) ` Rsn— Vi 1 w/( lo Rs’ Ron 


I 
i Rs+Ron R, n Vr RFR) 


= (Ren (Co +Ipn) * Rs +V1)/n > Vr ` (Rs +a») 


(9) 


3.2. Substitution of STC data into the I-V curve expressions 


Expressions (8) and (9) allow to carry out the calculations for 
determination of all the five unknown equivalent circuit para- 
meters, namely Ipn,Jo,Rs,Rsn, and n. The STC data is typically 
provided by the manufacturer and is utilized as follows. The short 
circuit current I$, the open circuit voltage Vš, and the maximum 
power point voltage-current pair (V, I},) of the STC I-V curve are 
referred to as basic points. The manufacturer usually supplies 
these data in a tabular form; moreover, these points can be easily 


extracted from the datasheet I-V curve with relatively high 
precision as well. 

Substitution of J, = 0 and V; = Vš. into (8) leads to the follow- 
ing open circuit STC relation: 

ak 
Voc = (U +p) ` Rsh—0° VF: w( b “Rep e+: ta - (10) 
VT 

Furthermore, substitution of V; = 0 and I, = /¥. into (9) leads to the 
subsequent short circuit STC expression, 


(+1) 
Rey ls ae 
r = Rolls py poy) OVE w -RelIRoy e i; a) 


Rs Rs V7 
The (Vý, ih) voltage-current pair determines the location of 
maximum power point on the STC I-V curve, 
= Vill): (12) 
Hence, combining (8) and (12), maximum power point STC 


relation is given by 


Ve =(54 Bn) Reh 


pk 
Th “(Rs Rsn) n° Vr y w( fo rE “Rh 5 ello In + Fon) Ron/n- a) i 
n Vr 
(13) 
The following condition determines the extremum of the power 
curve at I, =I: 


ddr: Vi) 


|. =" (14) 


hat, 


Combining of (8) and (14) leads to 


ABE Ron Zl (Re + Ray) + + W (Ratt fala vo) 
T 


T Rsh 
z z z 
n/N 7 Vi) ‘ ello In +55 Rsh /n Vr Ð) 


- n:-V;)}=0- (15) 
(am R i) 

To summarize, a system of four transcendental Eqs. (10), (11), (13), 

and (15) is now derived for STC. Note that the obtained set of 

equations is still insufficient for complete solution since there are 

five unknowns. 


3.3. Solving the system of equations 


A possible way to solve the derived system of equations is 
treating one of the unknowns as an independent parameter, which 
may be determined by fulfillment of an additional condition, 
namely minimization of deviation between the theoretical and 
measured l-V curves. Each one of the five unknowns may be 
selected as an independent one, though the ideality factor is 
probably the most convenient since its value resides in a known 
bounded range for any type of photovoltaic cells (1<7<2 for 
silicon cells/panels) [37]. 

Combining (1), (5) and (6) for open and short circuit STC 
conditions results in 


Bon Rsh — Vic 


k 
h= Ron < (eV3c/1: VÐ — 1) (16) 
And 

r ye Rsh +R; | Te(Rsh +Rs) — Vi. . eR /nVÐ —1 (17) 
phe sR Rh eVI — @URs/nV7)° 


respectively. Substitution of (16) and (17) into (13) and (15) results 
in the system of two transcendental equations with three 
unknowns Rs, Rsn and n, where ņ is an independent parameter, 
residing within a known range. The system of these two trans- 
cendental equations is solvable only by numerical methods. Each 
value of ņ creates a specific set of (Rs, Rsn), for which [5 and Bh may 
be found from (16) and (17). Each full set of five parameters 
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determines a specific J-V curve that passes through the three basic 
points and possesses a maximum power point at (Im, Vm) as shown 
in Fig. 2. Certain deviation exists between the curves except at the 
three basic points, where the curves must coincide. The correct 
value of 7 corresponds to the modeled curve with minimum 
deviation from the datasheet/measured I-V curve at STC. 

The minimization may be performed by visual comparison 
between modeled curves and the measured one as well as by the 
simplest error minimization algorithm. When a measured I-V 
curve is unavailable in raw data form and a datasheet picture only 
is available (which is a typical case), one can use digitizer software, 
e.g. [45], which helps digitizing scanned plots into a raw data with 
fair precision. 


3.4, Additional set of equations for thermal parameters calculation 


The first three steps allow a fairly accurate reproduction of the 
STC I-V curve. The final step is estimating I-V curves for arbitrary 
environmental conditions using temperature coefficients. The 
expression for the open circuit voltage temperature coefficient 
(TCVoc) can be obtained by combining the derivative of (8) with 
(2)-(4) and substituting T=T’, 


where 


ak 
Cw («® Ron) et amo (23) 
n 


T 


The values of E, and e can be then obtained by solving (20) 
and (22). 


4. Results and discussion 


The proposed method was first applied to Isofoton IS - 160 
mono-crystalline solar panel [46]. The seven parameters of the 
model (Rs, Rsn, IG, Bh n, Eg and e) were obtained in accordance 
with the four steps of the presented approach. Theoretical I-V 
curves derived using the proposed method for different environ- 
mental conditions were compared to experimental I-V data 
provided by manufacturer. All the equivalent circuit parameters 
were obtained from the manufacturer's data only. The results were 
obtained for the same environmental conditions as manufacturers’ 
and transferred into transparent I-V curve plots. These plots were 
superimposed with the manufacturers’ datasheets by correspond- 
ing scaling along each axis. 

The specifications of the IS - 160 panel and the manufacturer's 
I-V curve are shown in Table 1 and Fig. 3, respectively. Note that 


dV, EeloRsn + QRsn( 319 + IpnT enV nV7A(Eg — GRsn(lo 4 Bh 1,.)+ 3qnVi+qnV;7A) 
dT r-r nV3(A+ 1)qT* í 
where 

I Rsh pe pe r 
A = W Oeste Rn lo + En T,)/nV7) : (19) 

nV7 


The TCV,, is obtained from (18) by substituting I,=0 as 


TCV, EglRsh +4Rsh (GI + IpnT enV7 nV7B(Eg qRsn (lg t ph) t 3qn V7 t qnV7B) 
oc 


nql*V7(B+1) 
(20) 


where 


TOR sh Ra +1%,)/nVi 

B=W e sh (lo + pn)/" T) . (21) 
V7 

Similarly, the expression for the short circuit current temperature 

coefficient (TCI;-) may be derived from combining the derivative of 


(9) with (2)-(4) and substituting T=T’, 


(18) 


the data, provided by Isofoton in tabular form and the datasheet 
I-V curve slightly differ (common to different solar panels man- 
ufacturers). The deviation remains within the tolerance (+ 5%) 
declared by the manufacturer. 

Since the proposed method is very sensitive to input data, 
graphical data was used here rather than tabular one. The estima- 
tion procedure was carried out using Mathematica software [47]. 
The estimation results for STC are summarized in Table 2. 

The estimated parameters, given in Table 2, were used for 
creating modeled I-V and P-V curves for STC (Fig. 4) and different 
temperatures (Fig. 5). Good agreement may be observed, char- 
acterized by deviations of 0.1-0.5%. 

Furthermore, the proposed method was applied to three 
additional solar panels, specified in Tables 3-5. The estimation 
results are summarized in Tables 6-8. The corresponding I-V and 


_ di, _ Rs||Rsh : ( * + xo pk Tk xo pk Eg — 3nV7q—nV7qC Vr ) 
TClsc = |, oe VRTCI Elf +VGG +T" e) + q + h) RIR +RR E 22) 
A Table 1 
O, L, Specification of the IS - 160 solar panel [46]. 
& cei AA ee ey, p panel [46] 
f s Electrical characteristics, STC From datasheet From datasheet I-V 
4 / \ table curve 
L 
V 3 LZoom-A. p” À Number of series-connected 72 72 
ae cells ns 
8 Maximum power Pm, [W] 160 (+ 3%) 155.67 
R zB ——'- n-L3 Maximum power voltage Vm, 35 3532 
S Se eee IV] 
--- nel Maximum power current Im, [A] 4.57 4.42 
Z Shot circuit current Ise [A] 4.9 4.81 
Voltage (V7 Open circuit voltage Vo., [V] 43.8 43.2 
Temperature coefficient of Və», —0.378 —0.353 
10 40 30 mr 
Temperature coefficient of I,,, 0.0254 0.0254 
Fig. 2. The family of I-V curves with different ideality factors passing through the %IK 


basic points 
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Fig. 3. The STC I-V curve of IS - 160 panel, digitized from the datasheet. 


Table 2 
The estimated STC parameters for a single cell of the IS-160 panel. 


Parameter Estimated value 
Ideality factor, n 1.0 

Series resistance R,, [mQ] 8.0 

Shunt resistance R,p, [2] 3.25 
Photocurrent Ipn, [A] 4.824 
Reverse bias saturation current lo, [pA] 0.33 

Band gap energy Eg, [eV] 1.16 
Temperature coefficient of Ipn, [A/K] 254 


L-V Carve at standard conditions of measuring 
? 


Carrent, (A/ 
Power, (W) 


OF ERAN URRMMMRANDOUERRN DY 
Vottage, (V) 


Fig. 4. The I-V and power curves of IS-160 at STC. Simulation results (dotted) are 
added to the manufacturer's data (solid). 
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Fig. 5. The I-V and power curves of IS-160 for different temperatures. Simulation 


results (dotted) are added to the manufacturer's data (solid). 


Table 3 
Specification of the Suntech CNPV-280P solar panel [48]. 


Electrical characteristics, STC From datasheet I-V curve 


Number of series-connected cells ns 72 
maximum power Pm, [W] 280 
Maximum power voltage Vin, [V] 36.9 
Maximum power current Im, [A] 7.6 
Shot circuit current Ise [A] 8.2 
Open circuit voltage Voc, [V] 44.6 
Temperature coefficient of Voc, %/K -0.3 
Temperature coefficient of Ise %/K 0.05 
Table 4 


Specification of the Schuco SF-160-24-M175 solar panel [49]. 


Electrical characteristics, STC From datasheet I-V curve 


Number of series-connected cells n, 72 
Maximum power Pm, [W] 175 
Maximum power voltage Vm, [V] 36 
Maximum power current Im, [A] 4.86 
Shot circuit current Isc, [A] 5.2 
Open circuit voltage Voc, [V] 44.8 
Temperature coefficient of Vse %/K —0.48 
Temperature coefficient of Is, %/K 0.04 
Table 5 


Specification of the Kyocera KC200GT solar panel [50]. 


Electrical characteristics, STC From datasheet I-V curve 


Number of series-connected cells n, 54 
Maximum power Pm, [W] 200 
Maximum power voltage Vin, [V] 26.3 
Maximum power current Im, [A] 7.61 
Shot circuit current Ise [A] 8.21 
Open circuit voltage Voc, [V] 32.9 
Temperature coefficient of Vse %/K — 0.374 
Temperature coefficient of Ise %/K 0.039 
Table 6 


The estimated STC parameters for a single cell of the CNPV-280P panel. 


Parameter Estimated value 
Ideality factor, n 1.0 

Series resistance R,, [mQ] 4.2 

Shunt resistance Ren, [Q] 43 
Photocurrent Ipn, [A] 8.19 

Reverse bias saturation current lo, [nA] 0.27 

Band gap energy Eg, [eV] 1.19 
Temperature coefficient of Ipn, [A/K] 420 


power curves of these panels along with the manufacturers' 
supplied data are shown in Figs. 6-8. It can be concluded that the 
estimated curves nearly coincide with the manufacturer provided 
curves, indicating the feasibility and high precision of the proposed 
method. 


5. Conclusions 


An improved method for determining the parameters of single- 
diode photovoltaic cell equivalent circuit was developed in the 
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Table 7 
The estimated STC parameters for a single cell of the SF-160-24-M175 panel. 


Parameter Estimated value 
Ideality factor, n 1.2 

Series resistance R,, [mQ] 71 

Shunt resistance Ren, [2] 10.2 
Photocurrent Ipn, [A] 5.18 

Reverse bias saturation current Ig, [nA] 10.2 


Band gap energy Eg, [eV] 1.16 
Temperature coefficient of Ipn, [WA/K] 


Table 8 
The estimated STC parameters for a single cell of the KC200GT panel. 


Parameter Estimated value 
Ideality factor, n 1.05 

Series resistance R,, [mQ] 4.85 

Shunt resistance Ren, [2] 3.6 
Photocurrent Ipn, [A] 8.2 

Reverse bias saturation current lo, [nA] 1.05 

Band gap energy Eg, [eV] 1.22 
Temperature coefficient of Ipn, [A/K] 388 


I-V Curves of CNPY-280P at different cell temperature 
kW/m? 


IRRADIANCE: AML, 


Current(A) 


Voltage(V) 


presented work. The set of experimental data required for the 
extraction of model parameters included manufacturer provided 
STC data only. 

In order to calculate the five unknown parameters for STC 
conditions, the method combined a solution of algebraic equations 
system to obtain four unknown parameters; whereas the fifth 
parameter was obtained by minimization of the divergence 
between the modeled and experimental I-V curves. Two addi- 
tional “thermal” parameters were derived from the solution of 
another system of algebraic equations based on the solar cell/panel 
temperature coefficients. A number of I-V curve points besides the 
three basic points (short and open circuits, maximum power) were 
required to obtain an accurate estimation of the ideality factor. 
One can obtain these points directly from the STC I-V curve 
available in the cell/panel datasheet. For this purpose, dedicated 
digitizing software was applied. 

Since the accuracy of the model parameters strongly depends 
on the accuracy of the input data, the accuracy of the data offered 
by the manufacturer (typically 3%) is not always enough. The 
precision has to be in the 1-1.5% range for obtaining correct model 
parameters. High accuracy of the temperature coefficients is also 
necessary. Correction of manufacturer provided is often necessary 
and may be performed by deriving these parameters from the 
manufacturer provided graphical data. 


I-V Curves of CNPV-280P at different irradiance 
Cell Temp. 25°C 


9 
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Fig. 6. The I-V curves of CNPV-280P for different irradiance levels and temperatures. Simulation results (dotted) are added to the manufacturer's data (solid). 
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Fig. 7. The I-V curves of SF 160-24-M for different irradiance levels and temperatures. Simulation results (dotted) are added to the manufacturer's data (solid). 
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Fig. 8. The I-V curves of KC200GT for different irradiance levels and temperatures. Simulation results (dotted) are added to the manufacturer's data (solid). 


The proposed method was verified by comparison between 
modeled and datasheet I-V curves of several off the shelf panels of 
leading manufacturers. The results demonstrated average devia- 
tion of 0.1-0.5%, indicating the proposed method feasibility and 
high accuracy. 
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